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ABSTRACT 



Transport coefficients are obtained by incorporating a gauge principle into thermo 
field dynamics of inhomogeneous systems. In contrast to previous derivations, nei- 
ther imaginary time arguments nor perturbation theory in powers of a coupling 
. constant are used in the calculation. Numerical values arc calculated for the pion 

component in hot nuclear matter. 



1. Introduction 

In the presence of homogeneous matter, or in a heat bath, the irreducible rep- 
resentations of the space-time symmetry group are characterized by a continuous 
energy parameter rather than by a mass-shell constraint. It follows, that physical 
"particles" are unstable and therefore a naive perturbation theory starting from 
free on-shell particles may lead to inconsistencies. To circumvent this problem 
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£3 ; perturbative expansion in terms of generalized free fields may be defined, cf. ref. 

and papers quoted there. 

Thermo field dynamics (TFD, see ref. □) is a necessary input to this solution, 
since in contrast to the Schwinger-Keldysh (closed-time-path, CTP) method i of 
quantum statistical mechanics, it possesses two different (anti-) commuting repre- 
sentations of the underlying field algebra. Apart from this mathematical aspect 
however, TFD also contains a concept which extrexemly simplifies its practical 
application: Physical excitations are obtained as superposition of thermal quasi- 
particle states. Due to the facts outlined above, this must be a continuous super- 
position, with a weight function that is nothing but the spectral function of a field 
theory. 

To establish the notation, we first discuss TFD for a single bosonic oscillator, 
introducing canonical creation and annihilation operators a\ a, a', a for the two 
representations. The a, a operators annihilate the physical vacuum, and the two 
sets are transformed into another by means of an anti-unitary mapping, called the 
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tilde conjugation □. One may construct these representations in the Liouville space 
of the oscillator, where the generator of time evolution is 

H = ujN, N = a ] a-a ] a . (1) 

This operator possesses a continuous symmetry, since it is invariant under Bogoli- 
ubov transformations to a new set of operators 

where E is a (real) 2x2 matrix with determinant 1. These matrices form the 
symplectic group in two dimensions Sp(2) = SU(1,1), it is a noncompact and non- 
abelian group @. 

In thermo field dynamics, the new set of operators obeys thermal state condi- 
tions: 

£«W*))=0, £|W*)) = 0, (( ^||e tt = 0, ((W L \\? = 0, (3) 

where ((W || and HW^)) are left and right thermal state. It can be shown, that a 
particularly simple parameterization of TFD is obtained with 

= \\W)) , {W L \\ = ((1|| , B = ( 1 \ n ~ n ) . (4) 

Here, ||1)) is the Liouville vector associated with the unit operator, and \\W)) is 
the Liouville vector associated with the density matrix W of the quantum system. 
Statistical averages for an operator £ are calculated as ((1 cg)l || Vl^)) , i.e., as a matrix 
element instead of a trace over the Hilbert space of the oscillator. 

2. Gauging the Bogoliubov symmetry 

The Bogoliubov symmetry discussed above is a global symmetry, which for the 
simple case of a single oscillator means that it is time-independent. We now con- 
sider the question of a local Bogoliubov symmetry, which may arise in case the 
transformation matrix depends on some external parameters H. For convenience, 
we consider these external parameters as a vector r in some abstract space. 

Requiring the locality of a former global symmetry is called, in accordance with 
the methods applied to fundamental forces of nature, the gauging of the Bogoliubov 
symmetry. This name however is more than a simple analogy: It is known, that 
even in simple dynamical systems structures arise, that most closely resemble the 
coupling of the system to an external gauge potential i. 

The physical particle operators depend on the parameter vector r. Hence, the 
former will change with time, for infinitesimal changes of r according to 

r~r + S r: ( « « ) (1 + V{r)Sr) ( ) 



(_1^) T -(-l)VW- (5) 

V is a vector in r-space, P = (V 1 , V 2 , . . .), each component forming a 2 x 2-matrix. 

This implies, that the change of the physical particles with the external param- 
eters has a generator 

s-(_:i>)) T ^(^). (a) 

The total time evolution of the system is generated by the sum of H and <2, i.e., 
the Liouville equation for an operator £ reads 

i|l|£(t)»= (H+fi) ||£(t)». (7) 

Now consider the thermal quasi-particle operators, which under an r-evolution 
change with a different Q. From eqns. (@) and (|j) follows, that it is obtained 

as 

V\r) = B(r) V{r) B~ l {r) + B{r) (V r S _1 (r)) . (8) 

Hence, under a change of the operator description mediated by a Bogoliubov trans- 
formation, V transforms inhomogeneously like a gauge potential. The substitution 
of the ^-operators for the a-operators is therefore, in this new language, a gauge 
transformation. Finally, we require diagonality of the hamiltonian for thermal quasi- 
particles, i.e., 

(H + Q) [£] = u (e tt e - ?S) • (9) 

Returning to eqn. (||), this means that V has to vanish. This requirement fixes 
the connection between the physical particles and the statistical quasi-particles 
unambiguously: The change of the physical particle operators with the external 
parameters according to ([5]) has a gauge potential given by 

V(r) = - (V r B-\r))B(r) . (10) 

To summarize these considerations: Gauging the Bogoliubov symmetry of the time 
evolution in Liouville space, while keeping the thermal quasi-particle operators fixed, 
introduces a coupling of the physical particle representation to gradients in the 
external parameters. 

3. Interacting fields with local Bogoliubov symmetry 

An interacting charged scalar field at a fixed time can be expanded into momen- 
tum eigenmodes as 




and similarly for the tilde-field, k is the three-momentum of the modes, and in 
this notation a|_(t) creates a negatively charged excitation with momentum k, 
while a,k+{t) annihilates a positive charge. We will henceforth distinguish the two 
different charges by an additional index I = ± whenever possible. 

In general, the commutation relations of the a-operators are not known: They 
can be very complicated due to the interaction. However, for the purpose of cal- 
culating bilinear expectation values, e.g. propagators, it is sufficient to know the 
expectation value of the commutator functions. Hence, for this purpose, the scalar 
field operators may be considered as generalized free fields, with an expansion into 
modes with definite energy and momentum that are created and annihilated by a 
new set of operators. The details of the expansion procedure in general are outlined 
in refs. 0&i, its generalization to inhomogeneous states in 0>0. Here it is sufficient 
to write down the generalization of eqn. ([2]) 

a>ki(t) 

oo 

JdE jd 3 q (ME,(q + k)/2)) l/2 (Bf^E^k))* ( ) e" iEt . (12) 
o Eq 

Ai(E, k) is a positive weight function with support on the positive energy axis, and 
the transformation matrix has the form 

similar to the matrix in eqn. (|j). The Bogoliubov parameter iVj appearing here is 
the Fourier transform of a space-local quantity 



Nt(E, q, k) = JL Jd 3 z e-^'Q* m(E, (q + k)/2, 



(14) 



The thermal state conditions are similar to eqn. (|]) , but hold for both charges and 
for each value of E, k. The time derivative of the a-operators in the Heisenberg 
picture then is 

.d_ ( a k i(t) \ _ f a M (t) 
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which we use to determine the operator Q. The Z-factor is defined as 

oo 

Z kl = fdEAi{E,k), (16) 



and Qki — 1/(2^/). Obviously the additional term in the time evolution vanishes 
for translationally invariant n;, i.e., when Ni(E,k,q) is proportional to S 3 (k — q). 
This is consistent with eqn. (|T0|), henceforth the generator Q is 

It has a kernel 



4w 1 ' (17) 



^kl ^ql 

The dimension of the numerical factors in this kernel is mass 2 , due to the energy 
normalization factors contained in the a-operators for the interacting field (|TTD. 

4. Transport coefficients 

Observables are expressed as functionals of the interacting fields - and therefore 
as functionals of the operators a, afi. Hence, although the state we consider is 
stationary in terms of the basis defined by the ^-operators, momentum mixing as 
introduced above will result in a non-trivial time evolution of physical quantities. 
Two quantities we are interested in are the particle current operator of the complex 
scalar field and its expectation value, 



j(t, x) = i (0*V0 X - X V^) , j(t, x)=5( j(t, x) ) ( 19) 

and the energy current 

E(t, x) = (d^V^ + dtfaVfc) , E(t,x)=5{ E(t, x) ) . ( 20 ) 



The response of such a current to the gradients is, to first order in the gradients, 
given by 

t 

s(j(t,x)\=iJdr(\j(t-r,x),Q] \ . (21) 

t 



The commutator contains four field operators instead of two. The expectation 
value of this commutator therefore cannot be expressed completely through the 
commutator function of generalized free fields. Hence in the following, we neglect 
higher order correlations and use the canonical commutation relations of the £- 
operators to obtain in the stationary limit 

3l{t,x) = 

oo 

/fd 3 kc\ 3 k' 
dE J (2?r)3 Ai{E, fe) Ai{E, k') e i ( fc - fc )" (fc + k 1 ) Vi(k, k') (22) 

o 

for each species. The energy current is formally quite similar, and the total currents 
are for both cases 

j{t,x) = j + (t,x) - j_(t,x) 

E{t,x) = E+(t,x) + E_(t,x) . (23) 

In principle the above expressions can be calculated, when the spectral functions 
and the space-dependence of rii are given. However, it is better justified to perform 
a gradient expansion for the occupation number density n h since the usual spectral 
function is correct also only to first order in these gradients. The i-th. vector com- 
ponent of the /-charged currents generated by the inhomogeneity of the system is 
then 



j?(t,x) = 2«J-^^ r JdE (Ai(E,k)) x 

Lf'F' f dn l (E',k,x)dA l (E',k) \ 

j dEE \-^ ai&-f (24) 

, , f d 3 k k^ f 

E?(t,x) =^jj^y 3 ^r ( JdE E (Ai(E,k)) 2 x 



f dn l (E',k,x)dA l (E',k) \ 
dE \—d^) ■ (25) 

Before this is evaluated further, note that the expression in the curly brackets is 
nothing but the first term of a full gradient expansion of the product of ni and Au 
i.e., the Poisson bracket. It is therefore clear, that the above expression is useful 
also for situations, where the space-dependence of Ai cannot be neglected: the 
Poisson bracket then also has a contribution with a space-derivative acting on Ai 
and a momentum derivative acting on n\. The above expression is also interesting 
in view of the original equation (|2l| ), since the Poisson bracket is the analogon of 
the commutator, i.e., it contributes a factor h to the current. 

Using as local equilibrium spectral function the one for global equilibrium (cor- 
rect to first order in the gradients, cf. !), and as distribution function 

n±(E, k, x) = -kp; 3 *) = 1 ^ (26) 

±y } l-f ± (E,x) exp(P(x)(ETKx)))-l 



the derivative with respect to x has two components: Along the temperature gra- 
dient and along the gradient of \i. 

The same decomposition into T and /x-gradient terms then holds for the currents, 
i.e., with transport coefficients Lij we obtain 

VT 

j(t,x) =-LuV/i -L 12 — 

VT 

E(t,x) =-L 21 Vii -L 22 — . (27) 

Note, that L 12 and L 21 can be different, i.e. the Onsager relation Ly = Lji is not 
necessarily fulfilled. The reason for this is clear: in contrast to ordinary many- 
body quantum physics, the present formulation exhibits dissipation already on the 
tree-graph level S. In other words, the formulation of an interacting theory with 
continuous spectral functions is not micro-reversible, physical states have a finite 
lifetime. However, the complete system, including all species of particles and in- 
teractions, by definition has a single partition function. Hence, the global Onsager 
relations are fulfilled. 

From the above currents, one obtains the thermal conductivity for the interacting 
scalar field by subtracting the convective part 

Lu — d is the diffusion coefficient (or mobility). Before considering numerical 
results obtained with these equations, let us study the case of a simple spectral 
function 

ME, fc) - A- (E, k)^ {E2 _ ^ + ^ (29) 

with Vl\ = e\ + "fl and 7^ e^. In the absence of a chemical potential, L 12 and L 2 i 
are zero, but diffusion coefficient and thermal conductivity are, to lowest order in 
the width for \i = 0, 

I f d 3 k 1 d 
do = ~ W J^f^- k d-k tkn{tk){l + n{tk)) 

1 f d 3 k ktu d . . , . . . 

Here, n(efe) is the local Bose function taken at energy e^. Apart from the momentum 
factors, a similar representation for Ao was obtained in ref. 113. The difference can be 
attributed to the fact that in this reference a hydro dynamical picture was assumed 
together with <C 1. We thus make a comparison in the high temperature limit, 
where one obtains 

lim A = -T 3 . (31) 

T^oo 7 



Here, 7 is some momentum averaged width and A a numerical factor. The same 
limiting behavior is obtained in £3. 

The above expressions are also in accordance with naive expectations: The 
thermal conductivity and diffusion coefficient diverge, if the width of the particles 
is reduced, i.e., if the interaction is removed. In this case the time needed for the 
relaxation of a temperature disturbance is infinite. 

As stated above, The Poisson bracket contributes a factor h, while in a naive 
counting the inverse width contributes a factor 1/hto the transport coefficient. This 
cancellation seems to be the main source for the difficulties one has in obtaining 
kinetic coefficients from quantum field theory. 



5. Pionic transport properties 

In the next step, we apply the model to pions in a homogeneous gas of nucleons 
and A33-resonances. Instead of the above simplistic form, the more realistic pionic 
spectral function from ref. Q is used. It exhibits a strong mixing of the pion with 
a collective exciation of A-particle/nucleon-hole pairs, hence has a richer structure 
with two resonance branches. In the figure, the product AT and and the diffusion 
coefficient d are plotted for this spectral function. 

The high temperature behavior at the example density of p = 1.59x nuclear 
matter density inferred from the figure is 

A « 9897GeV/fm 2 x [T/GeV] 3 ' 70 , (32) 



i.e., it rises even faster than estimated by eqn. This is due to the decreasing 

width of the pion with temperature in the A-hole model - a fact, which was not 



included in eqn. (31). At lower temperatures, the thermal conductivity rises slower 
with temperature than the asymptotic expression. 

The density dependence of the coefficients A and d is quite small, with the general 
tendency to have a " stiffer" temperature dependence at lower baryon density. This 
is consistent with the fact, that the pion in the A-hole model becomes a free particle 
without the baryonic background. 

A comparison of these values to experimental data is beyond reach for the time 
being, we thus have to restrict the comparison to other calculations: The ther- 
mal conductivity in this field-theoretical calculation is substantially higher, than 
obtained in a fluid-dynamical picture fitted to pion scattering data at zero temper- 
ature and densityfcffil. 

This is consistent with the finding, that off-shell effects, which were treated 
consistently in the present approach, may substantially increase relaxation times El. 
Hence, within the framework of quantum field theory, the new method presented 
here gives reliable results for transport coefficients. 



Top panel: Thermal conductivity X temperature. 
Bottom panel: Diffusion coefficient. 

Thin straight lines: slopes OC T 3 , T 4 

Full thick line: A-hole model at p = 1.59 nuclear density; 
dashed line 0.92 and dash-dotted line 0.47 nuclear density. 

Dash-double-dotted in top panel 0, dotted line ref. El, 
both from 7T — 71 scattering data at T = 0, p = 0. 
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